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I. INTRODUCTION 



The neutrino induced reactions in nuclei at intermediate energies play an important role in the study of neutrino 
properties and their interaction with matter A good example of this is the search for neutrino oscillations, and 
hence physics beyond the standard model 0- Several experiments are planned or under construction 0, aimed at 
determining the neutrino oscillation parameters with high precision. The data analysis will be sensitive to sources 
of systematic errors, among them nuclear effects at intermediate energies (nuclear excitation energies ranging from 
about 100 MeV to 500 or 600 MeV), being then of special interest to come up with an unified Many-Body Framework 
(MBF) in which the electroweak interactions with nuclei could be systematically studied. Such a framework would 
necessarily include three different contributions: i) QE processes, ii) pion production and two body processes from 
the QE region to that beyond the A(1232) resonance peak, and iii) double pion production and higher nucleon 
resonance degrees of freedom induced processes. Any model aiming at describing the interaction of neutrinos with 
nuclei should be firstly tested against the existing data of interaction of real and virtual photons with nuclei. There 
exists an abundant literature on this subject, but the only model which has been successfully compared with data at 
intermediate energies and that systematically includes the first and the second of the contributions, and partially also 
the third one, mentioned above, is that developed in Refs. ^ (real photons), and 1] (virtual photons). This model 
is able to describe inclusive electron-nucleus scattering, total nuclear photo- absorption data, and also measurements 
of photo- and electro-nuclear production of pions, nuclcons, pairs of nucleons, pion-nucleon pairs, etc. The building 
blocks of this model are: 1) a gauge invariant model for the interaction of real and virtual photons with nucleons, 
mesons and nucleon resonances with parameters determined from the vacuum data, and 2) a microscopic treatment 
of nuclear effects, including long and short range nuclear correlations ,5i], FSI, explicit meson and A(1232) degrees 
of freedom, two and even three nucleon absorption channels, etc. The nuclear effects are computed starting from a 
Local Fermi Gas (LEG) picture of the nucleus, and their main features, expansion parameter and all sort of constants 
are completely fixed from previous hadron-nucleus studies (pionic atoms, elastic and inelastic pion-nucleus reactions, 
A— hypernuclei, etc.) The photon coupling constants are determined in the vacuum, and the model has no free 
parameters. The results presented in Refs. Q and are predictions deduced from the framework developed in 
Refs. 0-0 • One might think that LFG description of the nucleus is poor, and that a proper finite nuclei treatment 
is necessary. For inclusive processes and nuclear excitation energies of at least 100 MeV or higher, the findings of 
Refs. 0, Q and clearly contradict this conclusion. The reason is that in these circumstances one should sum 
up over several nuclear configurations, both in the discrete and in the continuum, and this inclusive sum is almost 
no sensitive to the details of the nuclear wave function, in sharp contrast to what happens in the case of exclusive 
processes where the final nucleus is left in a determined nuclear level. On the other hand, the LFG description of the 
nucleus allows for an accurate treatment of the dynamics of the elementary processes (interaction of photons with 
nucleons, nucleon resonances, and mesons, interaction between nucleons or between mesons and nucleons, etc.) which 
occur inside the nuclear medium. Within a finite nuclei scenarious, such a treatment becomes hard to implement, and 
often the dynamics is simplified in order to deal with more elaborated nuclear wave functions. This simplification of 
the dynamics cannot lead to a good description of nuclear inclusive electroweak processes at the intermediate energies 
of interest for future neutrino oscillation experiments. 

Our aim is to extend the nuclear inclusive electron scattering model of Ref. 0, including the axial CC degrees of 
freedom, to describe neutrino and antineutrino induced nuclear reactions. This is a long range project; in this work we 
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present our model for the QE region, and hence it constitutes the first step towards this end. We also present results 
for the inclusive muon capture in ^^C and predictions for the LSND measurement of the reaction ^^C {v^,n~)X near 
threshold. Both processes are clearly dominated by the QE contribution and are drastically affected by the inclusion 
of nuclear correlations of the RPA type. We find 
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= 3.2 X 10'' s"^ a{v^) = 11.9 X 10"''° cm^ (1) 



in a good agreement with data (discrepancies of the order of 10-15% for the muon capture rate), despite that those 
measurements involve extremely low nuclear excitation energies (smaller than 15-20 [25-30] MeV in the first [second] 
case), where the LFG picture of the nucleus might break down. However, it turns out that the present model provides 
one of the best existing combined description of these two low energy measurements, what increases our confidence on 
the QE predictions of the model at the hi^er transferred energies of interest for future neutrino experiments. Some 
preliminary results were presented in Ref. [3 • 

There exists an abundant literature both on the inclusive muon capture in nuclei, [8j-|l8|. and on the CC neutrino- 
nucleus cross section in the QE reg ion 0]-[33- Among all these works, we would like to highlight those included 
in Refs. ^3] (/^—capture) and |2lj ] (CC QE scattering) by Oset and collaborators. The framework presented in 
these works is quite similar to that employed here. Nicely and in a very simple manner, these works show the most 
important features of the strong nuclear renormalization effects affecting the nuclear weak responses in the QE region. 
The main differences with the work presented here concern to the RPA re-summation, which, and as consequence of 
the acquired experience in the inclusive electron scattering studies is here improved, by considering effects not 
only in the vector-isovector channel of the nucleon-nucleon interaction, but also in the scalar-isovector one. Besides, 
a more complete tensorial treatment of the RPA response function is also carried out in this work, leading all of these 
improvements to a better agreement to data. 

In addition here, we also evaluate the FSI effects for intermediate nuclear excitation energies, not taken into account 
in the works of Ref. |2l|. on the neutrino induced nuclear cross sections. 

The paper is organized as follows. After this introduction, we deduce the existing relation among the CC neutrino 
inclusive nuclear cross sections and the gauge boson W selfenergy inside the nuclear medium (Sect.m). In Sects. Hill 
and II VI we study in detail the QE contribution to the neutrino and antineutrino nuclear cross section, paying a special 
attention to the role played by the strong renormalization of the CC in the medium fSect. UlL^ and to the FSI effects 
('Sect. Ini The inclusive muon capture in nuclei and the relation of this process with inclusive neutrino induced 
reactions are examined in Sect0 Results and main conclusions of this work are compiled in Sects. IVII and IVIII 
Finally, in the Appendix, some detailed formulae are given. 

II. CC NEUTRINO INCLUSIVE NUCLEAR REACTIONS 

A. General Formulae 

We will focus on the inclusive nuclear reaction driven by the electroweak CC 

yi{k) + Az ^ I- {k') + X (2) 

though the generalization of the obtained expressions to antineutrino induced reactions, neutral current processes, or 
inclusive muon capture in nuclei is straightforward. 

The double differential cross section, with respect to the outgoing lepton kinematical variables, for the process of 
Eq. lO is given in the Laboratory (LAB) frame by 

dVL{k')dE[ |fc I 47r2 ^ 

with k and k' the LAB lepton momenta, E[ = {k'^ + mfy/^ and m; the energy, and the mass of the outgoing lepton 
(m^ = 105.65 MeV, rUe = 0.511 MeV ), G = 1.1664 x 10"" MeV"^^ ^he Fermi constant and L and W the leptonic 
and hadronic tensors, respectively. The leptonic tensor is given by (in our convention, we take £9123 — +1 and the 
metric g^'^' = {+,-,-,-)): 

Lf.a = + ii" - k'k, + k'^k^ - g^^k ■ k' + it^ao^pk^'k^ (4) 
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The hadronic tensor includes all sort of non-leptonic vertices and corresponds to the charged electroweak transitions 
of the target nucleus, i, to all possible final states. It is thus given by^ 

= ^E(2^)''^'(^}-^-9)(/bc^c(o)N)(/bc^c(o)N)* (5) 
* / 

with the four-momentum of the initial nucleus. Mi = P^ the target nucleus mass, Pj- the total four momentum 
of the hadronic state / and q = k — k' the four momentum transferred to the nucleus. The bar over the sum denotes 
the average over initial spins, and finally for the CC we take 

j^, = W„7''(l - 75)(cos0c*d + sm0c^s) (6) 

with and quark fields, and 9c the Cabibbo angle (cos^c — 0.974). By construction, the hadronic tensor 

accomplishes 

W^" = W^" + (J) 

with W^"^ (Wi"^) real symmetric (antisymmetric) tensors. To obtain Eq. Q we have neglected the four-momentum 
carried out by the intermediate M^— boson with respect to its mass, and have used the existing relation between the 
gauge weak coupling constant, g = e/ smOw, and the Fermi constant: = g^/8M^, with e the electron charge, 

9w the Weinberg angle and Mw the W^— boson mass. 

The hadronic tensor is completely determined by six independent, Lorentz scalar and real, structure functions 

— = -a^'''W^ + W9 + i- -^W^ + ^-^Wa + —^———^W,, + i— —^W„ (8) 

2M, 9 yvi+ "^'2 + ^ Mf * ^ 2Mi " ^ 2Mf ^ ^ ' 

Taking qin the z direction, ie, q = \q\uz, and P^ = (Mi,0), it is straightforward to find the six structure functions in 
terms of the M^™, W^^ = W^"^ ^ W^'^, W^^ and W'^^ components of the hadronic tensor^. After contracting with the 
leptonic tensor we obtain 



d^a^i \k'\E[M,G^ \ ^^^^ .^9' ^9' ^^E^ + E' 



' ' J 2Wi sin^ — + W2 cos^ W3 " ' sin^ - 



m 
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dn{k')dE'i 7r2 1 2 2 M, 2 E[{E[ + \k' 



Wi cos 9' cos 9' 

2 



^[+M _ cos U ^ K COS 9' + M(^[±M 3in^ 9']-W, 

2 M,; M, 2 \Mf ^ M? 2M,, 



(10) 



with El, the incoming neutrino energy and 9' the outgoing lepton scattering angle. The cross section does not depend 
on Mi, as can be seen from the relations of Eq. @, and also note that the structure function VFg does not contribute. 



B. Hadronic Tensor and the Gauge Boson Selfenergy in the Nuclear Medium 

In our MBF, the hadronic tensor is determined by the boson selfenergy, 11^ (g), in the nuclear medium. We 

follow here the formalism of Ref. 0, and we evaluate the selfenergy, SJ^(fc; p), of a neutrino, with four-momentum k 
and helicity r, moving in infinite nuclear matter of density p. Diagrammatically this is depicted in Fig. ^ and we get 

-iK{k;p) = / ^ur{k){-i^iiiD,^{q) (-in^^('?;p)) 'D,A<ir^ J[^l^T_^.^ (-i^) r,}uAk) (ii) 



^ Note that: (i) Eq. J3 holds with states normahzed so that {p\p') = {2n)'^2po5^ {p — p'), (ii) the sum over final states / includes an 

integration J ^2^)^25 ' each particle j making up the system /, as well as a sum over all spins involved. 
^ These relations read 

W^"" _ 1 f.,.oo Al')' n.... ...... o'?°o„„,oA _ -W^" 



Wi = - , W2 = Uy""-hVy''"= + ^^^(Vl/"-W^^)-2f--Re H/^M , W3 

M- 1 / \ Tm W^^ 

Wi = —^(W'' -W''), 1^5 = — Reiy''^-f-(iy"-iy^^) , We = — r-r- (9) 

2lgr |g| V l^l / \i\ 
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with D^a{q) — g^ia + liila / / (^q^ — M^,^ + ie^ , I\t^{q;p) is the virtual selfenergy in the medium, 

7l = 7''(1 ~ 75)7 ^'^d spinor normalization given by uu — 2m. Since right-handed neutrinos are sterile, only the 
left-handed neutrino selfenergy, Yi^{k]p), is not zero and obviously Y,^{k] p) — 'Y^^YT^ik; p). The sum over helicities 
leads to traces in the Dirac's space a thus we get 

^Ak:p)^^l4\,^^w^¥A- (12) 



V2M^J (2^)4 A:'2 - + ie 

The neutrino disappears from the elastic flux, by inducing one particle-one hole (Iplh), 2p2h ••• excitations, 
A(1232)— hole (Ah) excitations or creating pions, etc... at a rate given by 

T{k;p)^-^lmJ:4k;p) (13) 

We get the imaginary part of Si, by using the Cutkosky's rules. In this case we cut with a straight vertical line (see 
Fig. GJ the intermediate lepton state and those implied by the M^— boson polarization (shaded region). Those states 
are then placed on shell by taking the imaginary part of the propagator, selfenergy, etc. Thus, we obtain for fc" > 



ImS,(fc) = / -— lm{n^^iq;p)L,,} (14) 



with the Heaviside function. Since TdtdS provides a probability times a differential of area, which is a contri- 

bution to the {vi, I) cross section, we have 

da — r(k; p)dtdS = —-rrlmJ^^ik; p)dtdS — — —li'nSu(k;p)d'^r (15) 

\k\ 

and hence the nuclear cross section is given by 

a = — ^ / lmJ:^{k;p{r))d^r (16) 
\k\ J 

where we have substituted E,y as a function of the nuclear density at each point of the nucleus and integrate over the 
whole nuclear volume. Hence, we assume the Local Density Approximation, (LDA) which, as shown in Ref. Q, is an 
excellent approximation for volume processes like the one studied here. Coming back to Eq. H14|l we find 



dn{k')dk'° 




im(n^v^ + n;;,^) - l^^ Re(^^v^ - n^)}e(g°) (i?) 




FIG. 1: Diagrammatic representation of the neutrino selfenergy in nuclear matter. 
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and then by comparing to Eq. Q, the hadronic tensor reads 

Wr = -e(<z") (^^^ j ^ Im [IV>^ + If^] (g; p) (18) 

wr = -©(9°) (^^) / ["^^ ^ ('^^ ^) 

As we see, the basic object is the selfenergy of the Gauge Boson (W^^) inside of the nuclear medium. Following the 
lines of Ref . Q , we should perform a many body expansion, where the relevant gauge boson absorption modes would 
be systematically incorporated: absorption by one nuclcon, or a pair of nucleons or even three nucleon mechanisms, 
real and virtual meson (tt, p, ■ ■ ■) production, excitation of A of higher resonance degrees of freedom, etc. In addition, 
nuclear effects such as RPA or Short Range Correlations'^ (SRC) should also be taken into account. Some of the 
VF— absorption modes are depicted in Fig. [3 Up to this point the formalism is rather general and its applicability 
has not been restricted to the QE region. In this work we will focus on the QE contribution to the total cross section, 
and it will be analyzed in detail in the next section. 



W"^NN — ► NN 

W+ n ^ p W+ N — ► A, N* W"^N ^ N ji , Np, ... 





FIG. 2: Diagrammatic representation of some diagrams contributing to the W+— selfenergy. 



^ For that purpose we use an effective interaction of the Landau-Migdal type. 
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III. QE CONTRIBUTION TO n^;;:^(g;p) 

The virtual can be absorbed by one nucleon leading to the QE contribution of the nuclear response function. 
Such a contribution corresponds to a Iplh nuclear excitation (first of the diagrams depicted in Fig. [21). To evaluate 
this selfenergy, the free nucleon propagator in the medium is required. 

S{p ■,p) = {i> + M)G{p ; p), G{p ; p) = (^-^-^1^ + J^<5(/ - E{p))Q{kp - \p |)) (20) 

with the local Fermi momentum kpir) = (37r^p(r)/2)^/'^, M = 940 MeV the nucleon mass, and E{p) — \J + p ^ _ 
We will work on a non-symmetric nuclear matter with different Fermi sea levels for protons, A:^, than for neutrons, 
k"^ (equation above, but replacing p/2 by pp or p„, with p — pp + p„). On the other hand, for the W^pn vertex we 
take 

<p-p' ^p^q\j^M\n-J >^ u(p -A^)u{p) (21) 
with vector and axial nucleon currents given by 

^" = 2COS0C X (pYiq^)!" + '^^'v^^<T"''q}j , A" = C0s9cGA{q^) X (^7"75 + ^^^^q"l5^ (22) 

with m^r = 139.57 MeV. Partially conserved axial current and invariance under G-parity have been assumed to 
relate the pseudoscalar form factor to the axial one and to discard a term of the form {p^^ + p'^)75 in the axial 
sector, respectively. Invariance under time reversal guarantees that all form factors are real. Besides, Due to isospin 
symmetry, the vector form factors are related to the electromagnetic ones'* 

Fril') = I (Friq') - F-iq')) , pyFYiq') = \ {ppF^iq^) - PnF^{q^)) (24) 
and for the axial form-factor we use 

GA{q^) = Y^ 1^772^1, 5A = 1.257, A/^ = 1.049 GeV (25) 

(1 - q^/MXY 

With all of these ingredients is straightforward to evaluate the contribution to the W^"*"— selfenergy of the first diagram 
of Fig. El 

- \n^{q\ q) = - cos^ ec (^^) ' J ^.A^^iP^ q)Gip; Pn)G{p + q; Pp) (26) 
with the CC nucleon tensor given by 

A'"'(p, g) = Tr I (^2Fi^7'' - 2iAiy^a^"9„ - Ga {l^l, ~ ^^^l^^^s)) (l^ + ^ + 

X i^FYi^ -f 2iAiv ^^o^^qp - Ga (l""!^ + ^^^^q'l^^ ) (l^ + M) | (27) 

The Dirac's space traces above can be easily done and the nucleon tensor can be found in Sect. I A ll of the Appendix. 
Subtracting the divergent vacuum contribution in Eq. 126() . we finally get from Eqs. (|18|1 and (|19p 

X (-»r)i(5» + E(p)-B(p + <,-))/l"''(p,5)|,„,s,j, } (28) 



We use the parameterization of Galster and collaborators '32"| 



1 + T 1 + r ii„ IJ.,iT W-q^Mt, 



with T = -q'^/iM^, Md = 0.843 MeV, tip = 2.792847, ^J,„ = -1.913043 and A„ = 5.6. 
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The (fip integrations above can be analytically done and all of them are determined by the imaginary part of the 
relativistic isospin asymmetric Lindhard function, Uii{q, kp, kp). Explicit expressions can be found in Sect.^of the 
Appendix. 

Up to this point the treatment is fully relativistic and the four momentum transferred to the nucleus can be 
comparable or higher than the nucleon mass^. At low and intermediate energies, RPA effects become extreniely 
large, as shown for instance in Ref To account for RPA effects, we will use a nucleon-nucleon effective force |23| 
determined from calculations of nuclear electric and magnetic moments, transition probabilities and giant electric 
and magnetic multipole resonances using a non-relativistic nuclear dynamics scheme. This force, supplemented by 
nucleon-A(1232) and A(1232)-A(1232) interactions 0-0, was successfully used in the work of Ref. Q on inclusive 
nuclear electron scattering. In this latter reference a non-relativistic LFG is also employed. Thus, it is of interest 
to discuss also the hadronic tensor of Eq. (|28|1 in the context of a non-relativistic Fermi gas. This is easily done by 
replacing the factors M / E{p) and M/E{p + q) in Eq. H28|) by one. Explicit expressions can be now found in Sect. lUl 
of the Appendix. 

Pauli blocking, through the imaginary part of the Lindhard function, is the main nuclear effect included in the 
hadronic tensor of Eq. H28(l . In the next subsections, we will study different nuclear corrections to M^'"'. 

To finish this section, we devote a few words to the Low Density Theorem (LDT). At low nuclear densities the 
imaginary part of the relativistic isospin asymmetric Lindhard function can be approximated by 

ImURiq, k"F,kPp) « -npn^Jiq^ + M- E{q)) (29) 

and thus one readily finds^ 

^ui + Az^i-+x~Na^^+n^i-+p, N^A-Z (31) 

which accomplishes with the LDT. For future purposes we give in Sect. ^ of the Appendix the i^i + n + p 

differential cross section. 



A. RPA Nuclear Correlations 



When the electroweak interactions take place in nuclei, the strengths of electroweak couplings may change from their 
free nucleon values due to the presence of strongly interacting nucleons Indeed, since the nuclear experiments on 
P decay in the early seventies , the quenching of axial current is a well established phenomenon. We follow here the 
MBF of Ref. 0, and take into account the medium polarization effects in the Iplh contribution to the VF— selfenergy 
by substituting it by an RPA response as shown diagrammatically in Fig. |21 For that purpose we use an effective 
ph-ph interaction of the Landau-Migdal type 



V = Co|/o(p) + fl){p)TlT2 + 50(P)C?1CT2 + g'f){p)SlS2TlT2^ 



(32) 



where a and f are Pauli matrices acting on the nucleon spin and isospin spaces, respectively. Note that the above 
interaction is of contact type, and therefore in coordinate space one has V{ri,r2) oc 5{fi — r2)- As mentioned before, 
the coefficients were determined in Ref. 33!] from calculations of nuclear electric and magnetic moments, transition 
probabilities, and giant electric and magnetic multipole resonances. They are parameterized as 



fMr)) = ^ft 



p(o)- 



1 - 



PiO) 



(ex) 



where 



ft^ - 0.07 
ft^ - -2.15 



/o 



l{ex) 



./(in) 



0.45 
0.33 



(in) 

9o 



gf^ = 90 = 0.575 



fi 

g;(")=g;(^-)= 3^ = 0.725 



(33) 



(34) 



^ The only limitation on its size is given by possible quark effects, not included in the nucleon form-factors of Eqs. 1231 - 1251 . 
® The energy of the outgoing lepton is completely fixed once the Fermi distribution of the nucleons is neglected. Thus all structure 
functions Wi get the energy conservation Dirac's delta into their definition. Indeed, we have 

= ^L^^^-^Siq'' + M- E(q)) X A"'' 
8ME{q) ^ ^ " 



(30) 
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and Co = 380MeVfm^. In the S = 1 = T channel {aarf operator) we use an interaction with exphcit tt (longitudinal) 
and p (transverse) exchanges, which has been used for the renormalization of the pionic and pion related channels in 
different nuclear reactions at intermediate energies Thus we replace, 

3 

co9o{p)^i^2fif2 -> fifa J2 '^i^2^«r' '^^r = (Q^^Mi) + ^ mj)Vt{q)) (35) 
with q = q/\q \ and the strengths of the ph-ph interaction in the longitudinal and transverse channel are given by 

f2 f /A2-m2\^ ^2 
= ^{^p\ K2_ 2 2_ 2 + Cp^2, Kp^ 2500 MeV, rup = 770 MeV (36) 

\ p q J q ™p 

The SRC functions g[ and g[ have a smooth dependence [1,113, which we will not consider here'', and thus we 




FIG. 3: Set of irreducible diagrams responsible for the polarization (RPA) effects in the Iplh contribution to the selfenergy. 



9 



will take g'i{q) = g't{q) = g' = 0.63 as it was done in the study of inclusive nuclear electron scattering carried out in 
Ref. jjy, and also in some of the works of Ref. 6||. Note that, coffo ^-nd g'/^/mJ differ from each other in less than 
10%. 

We also include A(1232) degrees of freedom in the nuclear medium which, given the spin-isospin quantum numbers 
of the A resonance, only modify the vector-isovector {S = 1 = T) channel of the RPA response function. The ph-Ah 
and Ah-Ah effective interactions are obtained from Eqs. I|35|l and H36|l by replacing a S, f ^ T, where S, T are 
the spin, isospin iVA transition operators 5] and f f* = 2.13 /, for any A which replaces a nucleon. 

Thus, the V lines in Fig. O stand for the effective ph(Ah)-ph(Ah) interaction described so far. Given the isospin 
structure of the W^NN coupling, the isoscalar terms (/o and go) of the effective interaction can not contribute to the 
RPA response function. We should stress that this effective interaction is non-relativistic, and then for consistency 
we will neglect terms of order ©(p^/M^) when summing up the RPA series. 

To start with, let us examine how the axial vector term {GaJ^J5T+/2, where t± — iTy are the ladder isospin 
operators responsible for the n to p and p to n transitions, T~^\n >= 2\p >) of the CC axial current is renormalized. 
As mentioned above, we will only compute the higher density corrections, implicit in the RPA series, to the leading 
and next-to-leading orders in the p/M expansion. The nonrelativistic reduction of the axial vector term in the nucleon 
current reads 

GAUr'{p')^ri5Ur{p) = 2MGaxI (^5'^^* + 9^" "'^^^/^ + ' ' ') Y^r, i = 1, 2, 3 (37) 

with p' = p+q, and Xr a non-relativistic nucleon spin-isospin wave-function. In Eq. (|37l) there is a sum on the repeated 
index i and the dots stand for corrections^ of order O (p^/Af^,p'^/M^, In the impulse approximation, this 

current leads to a CC nucleon tensor^, 

A^'ip, qrj., vector = 8Af2 [AT + AT) , AT = G\g^^g'^^5^^, AT = -G\ {g^^g^^ ^ g^^g^^) ^^^±^1 (38) 

with i,j ~ 1,2,3 and there is again a sum for repeated indices. The tensor A'^'^{p,q)\^^f^^ vector '^^^ ^^^'^ obtained 
from the non-relativistic reduction of A^" (p, q) in Eq. ljAl|l . The Ai contribution comes from the leading operator 
— g'' W_|_/2, and involves the trace of G\a^a^ (Iplh excitation depicted in the first diagram of Fig.|3Jl. Let us consider 
first this simple operator and only forward propagating (direct term of the Lindhard function) ph-excitations. Taking 
into account the spin structure of this operator the scalar term /q of the effective interaction does not contribute 
either, and thus we are left with the spin-isospin channel of the effective interaction, V^'^ 'j\fj2Tif2- Let us now 
look at the irreducible diagrams consisting of the excitation of one and two ph states (first and second diagrams of 
Fig. |3J . The contribution of those diagrams to the W-selfenergy is 

« <p|Z±|n><n|^|p>|rr(aV-'") 

+ <p\'-±-\n>< n\T\p X p\f\n X n\^\p > ^ Tr{a'a')Tr{a''a^)ViV 

^ ^ k,l = l 

= Uiq, klXp) + 2C7(9, fc]^, A:?,)V:i) (39) 

The excitation of three ph states gives a contribution of U{2U)'^ Y^k ^ar^ar = U{2U)'^ {'li'lj^i' + (^y" - Qi'ijWt') to 
n^. Thus, the full sum of multiple ph excitation states, implicit in Fig. 13 leads to two independent geometric series, 
in the longitudinal and transverse channels, which are taken into account by the following substitution in the hadronic 
tensor (VF'"') 



6'^8M'^G\luaj{q,k''^,kl) 8M^G\lml^U{q,k'}.,kPp) 



q^q^ _ qiqj 



1 - 2C/(<z, fc^, fc^)Vz(g) 1 - 2{/(g, fc» , A;?,)l/*(g) 



.8M^GllmU(q, kl, fcP) ( _ „ + ^ ) (40) 

A ^'^^ \i_2U{q,k-,k^^)Vi{q)\^ |l-2[/(g,fc^,fc^)T4(g)P;^ ^ 



This is justified because taking into account the q— dependence leads to minor changes for low and intermediate energies and momenta, 

where this effective ph-ph interaction should be used. 

Note that q°/M is of the order Igp/M^. 

Keeping up to next-to-leading terms in the p/M expansion. 
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The factor 2 in the denominator above and that in Eq. H39(l comes from the isospin dependence, ti • r2, of the effective 
ph-ph interaction. Taking account of Ah and backward (crossed term of the Lindhard function) propagating ph 
excitations (see Fig.EJ, not accounted for by U is readily done by substituting 2U in the denominator by U{q, kp) — 
Un + C/a, the Lindhard function of Ref. 35], which for simpHcity we evaluate^" in an isospin symmetric nuclear 
medium of density p. The different couplings for N and A are incorporated in Un and J7a and then the same 
interaction strengths V/ and Vt are used for ph and Ah excitations Q . Taking q in the z direction, Eq. H40|l implies 
that the axial vector contribution to the transverse {xx,yy) and longitudinal {zz) components of the hadronic tensor 
get renormalized by different factors 1/|1 — U[q,kp)Vt{q)\^ versus 1/|1 — U{q, kF)Vi{q)\^. 

Let us pay now attention to the term A2 in Eq. H38() . it comes from the interference between the —g^^a'^T+ /2 and 
T+gr'^° ((T • {p+p')) /4M operators in Eq. 137|l . The consideration of the full RPA series leads now to the substitution^^ 



8M'^G\lmU{q, kl, fc^) ^ SM^^G^-i^^i^M^L^ (41) 
\q\ \q\ \^-U{q,kF)Vi[q)\^ 



in the Qz and zO components of the hadronic tensor W^'^ . 

Keeping track of the responsible operators, we have examined and renormalized all different contributions to the 
CC nucleon tensor A^'^ , by summing up the RPA series depicted in Fig. 13 The 00, Oz, zz, xx and xy components of 
the RPA renormalized CC nucleon tensor^^ can be found in Sect. I A 2l of the Appendix. As mentioned above, since 
the ph(Ah)-ph(Ah) effective interaction is non-relativistic, we have computed polarization effects only for the leading 
and next-to-leading terms in the p/M expansion. Thus, order O (kpp'^ /M'^ , kpf'^ /M'^ , kpq^/M) has been neglected 
in the formulae of the Sect. I A 2l of the Appendix. We have made an exception to the above rule, and since /iy could 
be relatively large, we have taken l^vFY\q\/M to be of order 0(0) in the p/M expansion. Finally, we should stress 
that the scalar-isovector term of the effective interaction (/') cannot produce Ah excitations and therefore, when this 
term is involved in the RPA renormalization, only the nucleon Lindhard function (Um) appears (see coefficient Cn 
in Eq. (|X9|) ). 

To finish this subsection we will discuss the differences between the medium polarization scheme presented here and 
that undertaken in Refs. . There is an obvious difference, since in these latter references the scalar-isovector term 

(/') of the ph-ph effective interaction was not taken into account. In addition, there are some differences concerning 
the tensorial treatment of the RPA response function. In the framework presented in this work we firstly evaluate the 
Iplh hadronic tensor and all sort of polarization (RPA) corrections to the different components of this tensor. In a 
second step we contract it with the leptonic tensor and obtain the differential cross section^^. The RPA corrections do 
not depend only on the different terms of the nucleon currents, but also on the particular component of the hadronic 
tensor which is being renormalized. Thus, as it is obvious, the RPA corrections, in general, are different for each of the 
terms ((F]^)^ {FYY^FYfY, G\, G^, GaGp, F^Ga and F1{Ga. with Gp = 2MGAl{rnl - q^) ) appearing in the CC 
nucleon tensor. Besides, for a fixed term, the polarization effects do also depend on the tensor component. Indeed, 
we have already mentioned this fact in the discussion of Eq. 140f) . where we saw that the axial vector contribution 
to the transverse {xx,yy) and the longitudinal (zz) components of the hadronic tensor get renormalized by different 
factors. 

In the works of Refs. 0, |^ the Iplh hadronic tensor, without polarization effects included, is first contracted 
with the leptonic one. This contraction is denoted as, up to global kinematical factors, I^P i'^ those references. 

In a second step, the authors of Refs. ^2.01 study the medium polarization corrections to |Tp. They find out 

different medium corrections for each of the terms of the CC nucleon tensor {{FYY' , • • • , F2 Ga)i as we do. However, 
for a fixed term, they cannot independently study the effect of the RPA re-summation in each of the different tensor 
components, since they are not dealing with the hadronic tensor itself, but with the contraction of it with the leptonic 
one. As a matter of example, to account for the RPA corrections to the axial vector-axial vector term, the following 
substitution is given in Refs. |il2i.i21i | 

^ {i\l-U{q,kp)Vt{q)\^ + i\l-U{q,kp)Vi{q)\^ ^^^^ 



The functions C/jv a-nd f^A ^^re defined in Eqs.(2.9) and (3.4) of Ref. l35l . respectively. Besides, note that in a symmetric nuclear medium 
Um = 217+ backward propagating ph excitation. For positive values of the backward propagating ph excitation has no imaginary 
part, and for QE kinematics (7a is also real. 

To evaluate the longitudinal contribution we use 2pz + qz = ("ip + q) ■ q/\q \ = (2£;(p) + g") /\q \ = 2Mq''/\q \ + 0{p'^ /M'^ ,q^ /M). 
Besides, the transverse part of the effective interaction does not contribute since {S^k ~ ^z<ik) (2p + q)k = 0- 
These are the needed components to compute the hadronic tensor PV^", when if is taken in the z direction. 

Note that the differential cross section is determined by the 00, Oz, zz, xx and xy components of W^" through their relation to the 
M^2=l,- -,5 structure functions. See Eqs. j5J and 1101 . 
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The above substitution can be recovered from Eq. (|40|l by contracting this latter equation with 5ij, and replacing 
2U U. Thus, Eq. is strictly correct, neglecting terms"^'' of order p/M, only for the contribution to jTp 
obtained from the contraction of the hadronic tensor with the g^i, term^^ of the leptonic one. The prescription of 
Eq. (|42|l is not correct for those contributions to ^ ^ |Tp arising from the contraction of the fc^fc<j + A^fc^ terms of 
the leptonic tensor with the axial vector-axial vector contribution of W^'^ . Note however that neglecting the bound 
muon three momentum^^ and up to terms of order p/M , Eq. H42|) is correct for the study of inclusive muon capture in 
nuclei, where it was first used by the authors of Refs. [lj,|2l|, and it is also reasonable for neutrino-nucleus reactions 
at low energies, where the RPA effects are more important. 



B. Correct Energy Balance and Coulomb Distortion Effects 

To ensure the correct energy balance in the reaction Q for finite nuclei, the energy conserving 5 function in Eq. H28|l 
has to be modified [l^l2H . The energies E{p) and E{p + q ) in the argument of the 5 function refer to the LEG of the 
nucleons in the initial and final nucleus. In the Fermi sea there is no energy gap for the transition from the occupied 
to the unoccupied states and hence ph excitations can be produced with a small energy, Q^^'^{r) — Ep{r) — Ep{r) . 
However, in actual nuclei there is a minimum excitation energy, Q — M{Az+i) ~ M{Az), needed for the transition 
to the ground state of the final nucleus. For instance , this Q value is 16.827 MeV for the transition ^^Cgs -^^^ Ngg 
and the consideration of this energy gap is essential to obtain reasonable cross sections for low-energy neutrinos. We 
have taken it into account by replacing 

q'^q'-{Q-Q'''"'{r)) (43) 

in the 5— function of the right hand side of Eq. (|28|l . 

The second effect which we want to address here is due to the fact that the charged lepton produced in the reaction 
of Eq. is moving in the Coulomb field of the nucleus described by a charge distribution pchir). In our scheme, 
we implement the corrections due to this effect following the semiclassical approximation used in Ref. pH . Thus, we 
include a selfenergy (Coulomb potential) in the intermediate lepton propagator of the neutrino selfenergy depicted in 
Fig. n We approximate this selfenergy inside the LEG by 

Ec = 2ft'°Vc(r), Vc{r) = -ATTa[^ dr'r'^p,h{r') + dr' r' p,h{r')^ (44) 

with a = 1/137.036 and the charge distribution, p^h, normalized to Z. The evaluation of the imaginary part of the 
u self-energy in the medium requires to put the intermediate lepton propagator on the mass shell. Following the 
Cutkosky's rules, and neglecting quadratic corrections in Vc, we find 



1 ^ -.Jik'^^El) 



ITT- 



'-Q{E[{r) - mi) (45) 



where E[ is the asymptotical outgoing lepton energy in regions where the Coulomb potential can be neglected, and 
the local outgoing lepton energy, E[{r), is defined by energy conservation 



E[ (r) + Vc (r) = ^Jmf+JC'^{r) + Vc (r) = E^ (46) 

Because of the Coulomb potential, the outgoing lepton three momentum, is not longer conserved, and it becomes 
a function of r, taking its asymptotical value, fc', at large distances. Therefore, q should also be replaced by a local 
function: q'{r) — k — JC'{r). Furthermore, from the d^k' integration in Eq. (|14|1 . and considering now the locality 
of the three momentum, we get from phase space a correction factor \K.' {r)\E[{r) /\k'\E[. This way of taking into 
account the Coulomb effects has clear resemblances with what is called "modified effective momentum approximation" 
in Ref. |36| | . The use of a plane- wave approximation in the interaction region is equivalent to the assumption that the 
Coulomb potential does not change the direction of the particles when they leave the nucleus. It should therefore not 



Medium renormalization effects are taken into account in ttiese terms by means of the substitution of Eq. I41i . 
Note that the axial vector— axial vector contribution to W^'^ is order 0{p^ /M^). 
This is an accurate approximation and we will also make use of it in Sect. |V] 
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FIG. 4: selfenergy diagram obtained from the first diagram depicted in Fig 

state in the ph excitation. 



I^Jby dressing up the nucleon propagator of the particle 



strongly alter an outgoing negatively charged lepton wave packet, which asymptotically is spherical, after it leaves 
the nucleus except by slowing it down and thereby changing the average radial wavelength and amplitude as the wave 
moves to larger r. As it is shown in i_36|, for total cross sections this procedure works very accurately for muons down 
to low energies. For low- energy electrons and positrons it is less accurate, and the use of the Fermi function F{Z, E[) 
(113) widely accepted in the literature. Anyway, Coulomb effects are small and they become relatively sizeable 
only for neutrino induced reactions near threshold and/or for heavy nuclei. 

To summarize the results of Subsections 1111 A1 and HIl Bl our final expression for the hadronic tensor is given by 



2A'P Jq \k'\E'i ^ J {2ttY E{p) E{p + q ') 

eik-Ar) ~ \P m\p + q' I - kUr))i-n)5 {q'" + E{p) - E{p + q')) A^^ip, qX^^EiP) } (47) 

with g'" = - (Q - g^^°(r)), q'{r) = k-IC'{r), and A^p^ given in Sect.Elof the Appendix. 



C. FSI Effects 



Once a ph excitation is produced by the virtual VF— boson, the outgoing nucleon can collide many times, thus 
inducing the emission of other nucleons. The result of it is a quenching of the QE peak respect to the simple ph 
excitation calculation and a spreading of the strength, or widening of the peak. The integrated strength over energies 
is not much affected though. A distorted wave approximation with an optical (complex) nucleon nucleus potential 
would remove all these events. However, if we want to evaluate the inclusive {vi, l~) cross section these events should 
be kept and one must sum over all open final state channels. 

In our MBF we will account for the FSI by using nucleon propagators properly dressed with a realistic selfenergy 
in the medium, which depends explicitly on the energy and the momentum |39l |. This selfenergy leads to nucleon 
spectral functions in good agreement with accurate microscopic approaches like the ones of Refs. 0j The 
selfenergy of Ref. [s^ has a proper energy-momentum dependence plus an imaginary part from the coupling to the 
2p2h components, which is equivalent to the use of correlated wave functions, evaluated from realistic NN forces and 
incorporating the effects of the nucleon force in the nucleon pairs. Thus, we consider the many body diagram depicted 
in Fig. ^ (there the dashed lines stand for an NN interaction inside of the nuclear medium _5, 39]). However, a word 
of caution is in order since the imaginary part of this diagram presents a divergency. The reason is that when placing 
the 2p2h excitation on the mass shell through Cutkosky rules, we still have the square of the nucleon propagator with 
momentum p + q in the figure. This propagator can be placed on shell for virtual PF— bosons and we get a divergence. 
The divergence is not spurious, in the sense that its meaning is the probability per unit time of absorbing a virtual 
by one nucleon times the probability of collision of the final nucleon with other nucleons in the infinite Fermi sea 
in the lifetime of this nucleon. Since this nucleon is real, its lifetime is infinite and thus the probability is infinite, as 
well. The problem is physically solved [42| by recalling that the nucleon in the Fermi sea has a selfenergy with an 
imaginary part which gives it a finite lifetime (for collisions). This is taken into account by iterating, in the Dyson 
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equation sense, the nucleon selfenergy insertion of Fig. [S] in the nucleon hne, hence substituting the particle nucleon 
propagator, G(p; p), in Eq. H2UI) by a renormahzcd nucleon propagator, GpsiiP'i p)i including the nucleon selfenergy in 
the medium, Y,(jP ,p; p), 



Gfsi(p; p) 



-E{p)-^{p^,p-p) 



(48) 



with E[p) = M + p^/2M. As mentioned above, we use here the nucleon selfenergy model developed in Ref. l39ll . 
which led to excellent results in the study of inclusive electron scattering from nuclei jl| . Since the model of Ref. (23 
is not Lorentz relativistic and it also considers an isospin symmetric nuclear medium, we will only discuss the FSI 
effects for nuclei with approximately equal number of protons and neutrons, and using non-relativistic kinematics for 
the nucleons (see Sect. Oof the Appendix). Thus, we have obtained Eq. H48|l from the non-relativistic reduction of 
G(p; p), in Eq. (|2DJ), by including the nucleon selfenergy. 

Alternatively to the nucleon selfenergy language, one can use the spectral function representation 



Gfsi(p;p) = 



p^ ^ Lu — ie 



dco 



Sp{uj,p;p) 



uj + ie 



duj 



where Sh , Sp are the hole and particle spectral functions related to nucleon selfenergy S by means of 



Sp,h{^,p;p) 



lmT,{uj,p-p) 



Lj-E{p)~ReJ:{u,p;p) 



lmJ:{uj,p; p) 



with Lo > p or iv < fj, for Sp and Sh^ respectively. The chemical potential p is determined by 

By means of Eq. I|49|) we can write the ph propagator or new Lindhard function incorporating the effects of the 
nucleon selfenergy in the medium, and we have for its imaginary part (for positive values of q") 



(49) 



(50) 



(51) 



9(g°) 
47r2 



d^p 



dujSh{uj,p; p)Sp{q° + uj,p + q:p) 



(52) 



Comparing the above expression with that of the ordinary imaginary part of the non-relativistic Lindhard function, 
Eq. I|C3I) . one realizes that to account for FSI effects in an isospin symmetric nuclear medium of density p we should 
make the following substitution 



26(9°) 



d^p 
(2^ 



e{kUr) - \P m\p + q I - fc^(r))(-^)<5(g" + E{p) - E{p + q M'-'^ip, q)\pO=Eip) 



47r2 



dp dujSh{uJ,p;p)Sp{q + uj,p + cf; p)A''''{p,q)\pO^ 



:E{p) 



(53) 




FIG. 5; Insertion of the nucleon selfenergy on the nucleon line of the particle state. 
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in the expression of the hadronic tensor (Eq. pSfl l. The (Pp integrations have to be done numerically. Indeed, the 
integrations are not trivial from the computational point of view, since in some regions the spectral functions behave 
like delta functions. We use the spectral functions calculated in Ref. [s^j but since the imaginary part of the nucleon 
selfenergy for the hole states is much smaller than that of the particle states at intermediate nuclear excitation energies, 
we make the approximation of setting to zero ImS for the hole states. This was found to be a good approximation in 
|43| . Thus, we take 

Shito,p; p) = 5{lo - E{p))Q{ii - E{p)) (54) 

where Eiji) is the energy associated to a momentum p obtained self consistently by means of the equation 

E{p) = E{p) + Re^{E{p),p- p) (55) 

It must be stressed that it is important to keep the real part of S in the hole states when renormalizing the particle 
states because there are terms in the nucleon selfenergy largely independent of the momentum and which cancel in 
the ph propagator, where the two selfenergies subtract. 

On top of the FSI corrections examined here, one should also take into account the nuclear corrections studied 
previously in Subsections IIII Al and IIII Bl 

IV. CC ANTINEUTRINO INDUCED NUCLEAR REACTIONS 

The cross section for the antineutrino induced nuclear reaction 

vi{k)+ Az^l+{k')+X (56) 
is easily obtained from the expressions given in Sects. ITll and lllll with the foUowings modifications: 

• Changing the sign of the parity-violating terms, proportional to VF3, in the differential cross section, Eq. (|10|l . 

• Replacing the M^+selfenergy in the medium, H^jJ^, by that of the boson (H^). This is achieved by exchanging 
the role of protons and neutrons in all formulae, 11^ {pp{r), Pnir)) — {pn{r), pp{r)). 

• Changing the sign of Vc, which turns out to be repulsive for positive charged outgoing leptons. 

LFG 

• Correcting the LFG energy balance with the difference Q — Q (r), with Q — M{Az-i) — M{Az) and 
Q''''''{r) = E^{r)-EUr). 

V. INCLUSIVE MUON CAPTURE IN NUCLEI 

In this section we study the /i— atom inclusive decay, it is to say the reaction 

i^^-^'Xln,^Mk)+X (57) 

It is obvious that the dynamics that governs this process is related to that of antineutrino (Eq. ]b<i\i ) and neutrino 
(Eq. |j2Jl) induced nuclear processes, but a distinctive feature is that the nuclear excitation energies involved in the 
/i— atom decay are extremely low (smaller than w 20 MeV). In this energy regime one might expect important un- 
accuracies in the LFG description of the nucleus. However and due to the inclusive character of the process, we will 
see that our MBF leads to reasonable results, with discrepancies of the order of 10-15% at most, with RPA effects as 
large as a factor of two. We should emphasize that similar conclusions were achieved in the works of Ref. which 
also use a LFG picture of the nucleus. 

The evaluation of the decay width for finite nuclei proceeds in two steps. In the first one we evaluate the spin 
averaged decay width for a muon at rest^^ in a Fermi sea of protons and neutrons with N Z. In this first step, 
the strong renormalization effects (RPA) will be also taken into account and thus we will end up with a decay width. 



^'^ In what follows, we will neglect the three momentum of the bound muon. 
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f which will be a function of the proton and neutron densities. In the second step, we use the LDA to go to finite 
nuclei and evaluate 

''dV|0i,(f)|2f(pp(r),p„(r)) (58) 

where (t>is{T^) is the muon wave function in the Is state from where the capture takes place. It has been obtained 
by solving the Schrodinger equation with a Coulomb interaction taking account of the finite size of the nucleus and 
vacuum polarization |38|. Equation H58|l amounts to saying that every bit of the muon, given by the probability 
\'Pis{r)\'^d^r, is surrounded by a Fermi sea of densities pp{r), pn{r). The LDA assumes a zero range of the interaction, 
or equivalently no dependence on q. The q dependence of the interaction is extremely weak for the /i— atom decay 
process and, thus the LDA prescription becomes highly accurate [T^. 

The spin averaged muon decay width, in an infinite nuclear matter of densities Pp{r) and Pn{T), is related to the 
imaginary part of the selfenergy (see Fig. EJ, (Pp(^)j Pn(f )), of a muon at rest and spin r in the medium by 

r(Pp(0>P"(r-)) = —ImT.^ {pp{r),pn{r)) , = ^ XI (^^) 

^ r 

To evaluate the imaginary part of the selfenergy associated to the diagram of Fig. El the intermediate states are 
placed on shell in the integration over the internal variables. These states are those crossed by the dash-dotted line 
in Fig. El The evaluation of the /u— selfenergy is almost identical to that of a neutrino in Fig. ^ and thus we obtain^^ 
from Eq. lfTI|l . 

f . --L_|L. I.,{n- (,;.,W,„.W)L„} (60) 

with q'^ = — \k \ and \q \ = |fc |. For kinematical reasons, only the QE part of the selfenergy will contribute 

to the muon decay width and thus we find 

r(Pp(r),p„(r)) = / — — ^— T-L^,,T'''' (g;pp,p„) 



J (27r)3 2|fc| 

G^cos^^c /■+°°r2/ , , h , _^ "^L _^ , \ .r, 

\k\h + -^t4, + mp,tz \ d\k\ (61) 



27r^ L I " 2 ' ' 2 




There is a factor 1/2 of difference coming from the averaged over initial spins of the muon, besides the W —selfenergy arises since the 
negative muon decay process is related to the antineutrino induced process, and the contribution of parity violating terms flips sign. 
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Nucleus Rp [fm] 


Rn [fm] 


a [imj 


Q [MeV] Q [MeV] 


[MeV] 




1.692 


1.692 


1.082 


16.827 


13.880 


0.100 


16(3 


1.833 


1.833 


1.544 


14.906 


10.931 


0.178 


18q 


1.881 


1.975 


1.544 


1.144 


14.413 


0.178 




2.773 


2.81 


0.54 


3.546 


4.887 


0.336 


*«Ca 


3.51 


3.43 


0.563 


13.809 


1.822 


1.064 


**Ca 


3.573 


3.714 


0.563 


3.142 


6.170 


1.063 


^5 As 


4.492 


4.64 


0.58 


0.353 


1.688 


2.624 


"^Cd 


5.38 


5.58 


0.58 


2.075 


4.462 


4.861 


208 p^, 


6.624 


6.890 


0.549 


2.368 


5.512 


10.510 



(*) The parameter a is dimensionless for the MHO density form. 



TABLE I: Charge {Rp,a), neutron matter {R„,a) density parameters, Q,Q— values and negative muon binding energies for different 
nuclei. For carbon and oxygen we use a modified harmonic oscillator (MHO) density, p(r) = po(l + exp(— (r/ij)'^), while for the 

rest of the nuclei, a two-parameter Fermi distribution, p(r) = po/(l + exp((r — R)/a)), was used. 

where the tensor T^^^ is defined as 

X (-7r)5(qO + E{p) - E(p + q ))A'^WiP^ q)\pO=Em } 

= ti^^-^ + ta/"^" + it3e'"'"''/ag/3 + t4g^(Z" + i5 (^^g" + ^"g^) , with = (1,0). (62) 

The similitude of the above equation with Eq. (|28() is clear. As in this latter case, the cPp integrations in Eq. H62() 
can be done analytically (see Sect. ^ of the Appendix) and all of them are determined by the imaginary part of the 
relativistic isospin asymmetric Lindhard function, UrIq, kp, kp). For a non-relativistic Fermi gas, the decay width is 
easily obtained from Eq. (|62|l by replacing the factors M / E{p) and AI/E{p + q) by one. Analytical expressions can be 
now found in the Sect.[niof the Appendix. FSI effects can be also taken into account by performing the substitution 
of Eq. 

Thus, both the muon decay process in the medium and the electroweak inclusive nuclear reactions i'i{k) + Az — *■ 
l~{k') + X in the QE regime are sensitive to the same physical features, VF^pn— vertex, and RPA and FSI effects. 
However, in the muon-atom decay only very small nuclear excitation energies are explored, 0-25 MeV, while in the 
latter processes higher nuclear excitation energies can be tested by varying the incoming neutrino momentum. 

The Is muon binding energy, B^^ > 0, can be taken into account, by replacing — > = — Bjf . This 
replacement leads to extremely small (significant) changes for light (heavy) nuclei, where Bjf is of the order of 0.1 
MeV (10 MeV), see Table U 

Finally, the correct energy balance in the decay can be enforced in the LFG by replacing 

9° ^ g° - (Q - Q^^'^ir)) - - jfc I - (Q - Q'^'^'^ir)) (63) 

in Eqs. l|HT|) and 

VI. RESULTS 

Firstly, we compile in Table HI the input used for the different nuclei studied in this work. Nuclear masses and 
charge densities are taken from Refs. (4^] and respectively. For each nucleus, we take the neutron matter density 
approximately equal (but normalized to N) to the charge density, though we consider small changes, inspired by 
Hartree-Fock calculations with the density-matrix expansion '46'| and corroborated by pionic atom data l47l| . However 
charge (neutron) matter densities do not correspond to proton (neutron) point-like densities because of the finite 
size of the nucleon. This is taken into account by following the procedure outlined in Section 2 of Ref. (see Eqs. 
(12)-(14) of this reference). 
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LDT 


Pauli+Q 


RPA 


SM \15] 


SM [27] 


CRPA [IS] 




Exp 
















LSND'95 [50] 


LSND'97 [51] 


LSND'02 [52] 


a 66.1 


20.7 


11.9 


13.2 


15.2 


19.2 


8.3 ± 0.7 ± 1.6 


11.2 ± 0.3 ± 1.8 


10.6 ± 0.3 ± 1.8 














KARMEN [53] 


LSND [54] 


LAMPF [55] 


a {ue,e~) 5.97 


0.19 


0.14 


0.12 


0.16 


0.15 


0.15 ±0.01 ±0.01 


0.15 ±0.01 ±0.01 


0.141 ±0.023 



TABLE II; Experimental and theoretical flux averaged ^"^Ciy^.ji )X and ^^C{ug,e )X cross sections in 10 cm^ units. We label 
our predictions as in Fig. |7| We also quote results from other calculations (see text for details). 



A. Inclusive Neutrino Reactions at Low Energies 



In this subsection we present results obtained by using non-relativistic kinematics for the nuclcons. We do not 
include FSI effects, since in Sect. IIII O we made the approximation of setting to zero ImE for the hole states. For low 
nuclear excitation energies (< 60), this approximation is not justified, because the imaginary part of the selfenergy 
of particle and hole states are comparable [s^. The inclusion of FSI effects would lead to a quenching of the QE 
peak of the bare ph calculation, and a spreading of the strength [1I I2M lisL l49l | . However FSI effects on integrated 
quantities are small. From the results of the next subsection we estimate in ~ 5-10% the theoretical error of the 
integrated cross sections and total muon capture rates presented in this subsection. 

The processes studied in this subsection explore quite low nuclear excitation energies (< 25-30 MeV), and hence 
one might expect that a proper finite nuclei treatment could be in order. Indeed, these processes are sensitive to 
the excitation of giant resonances 

El III El 113 ■ As mentioned in the Introduction, our purpose is to describe the 
interaction of neutrinos and antineutrinos with nuclei at higher energies (nuclear excitation energies of the order of 
100-600 MeV) of interest for future neutrino oscillation experiments. However, our model provides a good description 
of the low energy inclusive measurements analyzed in this subsection. RPA correlations play an essential role and 
lead to reductions as large as a factor of two. We should remind that the effective interaction appearing in the RPA 
series was fitted in Ref. |83j to giant resonances, and thus our approach incorporates the mechanism which produces 
those resonances in finite nuclei. 

At low energies, finite nuclei effects are expected to be sizeable for outgoing lepton energy distributions. There exist 
discrete and resonance state peaks, and the continuum distribution significantly differs from the LFG one. However, 
the integrated strength over energies, including the discrete state and resonance contributions, remains practically 
unchanged, which explains the success of our model to describe integrated-inclusive magnitudes. A clear example 
of this can be found in Ref. "4^ where the inclusive decay width of muonic atoms by using a shell model with final 
neutron states lying both in the continuum and in the discrete spectrum are calculated. The results are compared 
with those obtained from a LFG model. Both models^^ are in quite good agreement within a few percent when the 
shell model density is used in the LFG calculation. Being an integrated, inclusive observable, the total capture width 
is quite independent of the fine details of the nuclear wave functions. Similar conclusions were reached in the study 
of the radiative pion capture in nuclei, {Az — ''r~)i5oy„(j 7 ± AT, performed in Ref. There, the predictions of 

a continuum shell model were also extensively compared to those deduced from a LFG picture of the nucleus. The 
differences found, among the integrated decay widths predicted by both approaches, were, at most, of the order of 
4% (see Table 5 of first entry in Ref. 49]). 



1. The Inclusive Reactions C(z/^,/i )X and C{ue,e )X Near Threshold 



In order to compare with the experimental measurements we calculate flux averaged cross sections 



j \ dEMEu)WiE,), j [ W{E,)dE,, (64) 



In the LSND experiment at Los Alamos, the inclusive ^^C{v^, ix~)X cross section was measured using a pion decay 
in flight beam, with energies ranging from zero to 300 MeV, and a large liquid scillantor detector |50l Isil Is^ . The 
muon neutrino spectrum, W{E^), is taken from Ref. ^Ed\ and it is plotted in the left bottom panel of Fig. {7\ We 



For simplicity in the calculations of Ref. i48i . RPA effects are not considered and the static form of the nucleon CC current is employed. 
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FIG. 7: Predictions for the LSND measurement of the ^^C {Uf^, fi~ )X reaction (left panels) and the ^-^C {fe, e~)X reaction near threshold 
(right panels). Results have been obtained by using non-relativistic kinematics for the nucleons and without FSI. Top: Vf^ and cross 
sections multiplied by the neutrino fluxes, as a function of the neutrino energy. In addition to the RPA calculation (solid line), we show 
results without RPA correlations and Coulomb corrections (dotted line), and also (dashed line) the low density limit of Eq. Middle: 
Differential muon and electron neutrino cr oss s ections at Ei,^ = 179.5 MeV (left) and i?^^ = 46.2 MeV (right), as a function of the energy 
transfer. Bottom: Neutrino spectra from l5d| (left) and Eq. I65i (right). 
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fix i?™'" and Ef^'^^ to 123.1 and 300 MeV, respectively. The electron neutrino beams used in experiments (LAMPF, 
KARMEN, etc.) have relatively low energies. Such neutrinos do not constitute a monochromatic beam, and their 
spectrum^" is plotted in the right bottom panel of Fig. |3 The bare ph strength spreading due to the FSI might affect 
the inclusive, flux averaged cross-section because of the energy variations in the neutrino flux. As an illustration, if 
some of the strength is shifted to higher energies then some of the low-energy neutrinos will not be able to excite 
it, compared with the case when the strength is not spread out. Of course these effects are not very large, because 
some strength is also moved to lower energies and compensates this. These uncertainties contribute to the ~ 5-10% 
theoretical error mentioned above. 

Our results for the ^^C{i'f^, ^~)X and ^^C(i^e,e~)X reactions near threshold are presented in Fig. [7| and Table ITU 
As can be seen in the table, the agreement to data is remarkable. Nuclear effects turn out to be essential, and thus 
the simple prescription of multiplying by a factor of six (the number of neutrons of ^^C) the free space vin —> pl~ 
cross section overestimates the flux averaged cross sections by a factor of 5 and of 40 for the muon and electron 
neutrino induced reactions, respectively. The inclusion of Pauli blocking and the use of the correct energy balance 
in the reaction lead to much better results, but the cross sections are still badly overestimated. Only once RPA and 
Coulomb corrections are included a good description of data is achieved. RPA correlations reduce the flux averaged 
cross sections by about a factor of two, while Coulomb distortion significantly enhances them, in particular for the 
electron neutrino reaction where this enhancement is of about 30%. 

In Table m a few selected theoretical calculations (large basis shell model (SM) results of Refs. [Tsl and the 
continuum RPA (CRPA) ones from Ref. ,1^]) are also quoted. Our approach might look simplified with respect to 
the ones just mentioned, but in fact it is also an RPA approach built up from single particle states of an uncorrelated 
Fermi sea. This method in practice is a very accurate tool when the excitation energy is sufficiently large such 
that relatively many states contribute to the process. Obviously, because of its nature, the method only applies to 
inclusive processes and it is not meant to evaluate transitions to discrete states. The adaptation of the method to 
finite nuclei via the LDA has proved to be a rather precise technique to deal with inclusive photonuclear reactions 
and response functions in electron scattering The effective ph(Ah)-ph(Ah) interaction used in the RPA series has 
been successfully employed in different processes P, 0, . There are two distinctive features of this interaction in the 
S = T — I channel, which are not incorporated in most of the finite nuclei approaches: i) it incorporates explicit 
pion and rho exchanges and thus the force in this channel is splited into longitudinal and transverse parts, and ii) 
it includes resonance A degrees of freedom. The inclusion of Ah components in the RPA series reduces the LSND 
flux averaged ^^C {v^t^~)X cross section by about a 15%, while the reduction factor is about four times smaller for 
the electron neutrino reaction, because in this latter case, the larger contributions to the flux averaged ^^C(z^e, e~)X 
cross section comes from very low (< 20 MeV) nuclear excitation energies (see Fig.[7)|. In addition, a correct tensorial 
treatment of the RPA hadronic tensor is also important, and it explains the bulk of the existing differences between 
our results and those obtained in Ref. (see Subsect. IIII Al for details). As a matter of example, in Ref. |23| a 
value of 16.7 ± 1.4 x 10^**° cm^ is predicted for the LSND flux averaged ^^C(f^,/i~)X cross section. This value is 
about a 40% higher than our result, despite of using quite similar ph(Ah)-ph(Ah) effective interactions. Differences 
are significantly smaller for the electron neutrino flux averaged cross section, since this reaction is sensitive to quite 
lower energies. 

In the middle panels of Fig.[7| we plot the outgoing lepton energy distribution for an incoming neutrino energy near 
the maximum of (t{E^)W{E^) (top panels). We see in these plots the range of energies transferred [E^ — E[ — Q) 
to the daughter nucleus: 25-30 MeV for the muon neutrino reaction and less than 10 MeV for the electron neutrino 
process. Finite nuclei distributions will present some discrete state and narrow resonance peaks, but the integrated 
strength over energies would not be much affected though, as we have already discussed. 

2. Total Nuclear Capture Rates for Negative Muons 

After the success in describing the LSND measurement of the reaction ^^C {v^, l^~)X near threshold, it seems natural 
to further test our model by studying the closely related process of inclusive muon capture in ^^C. Furthermore, and 
since there are abundant and accurate measurements of nuclear inclusive muon capture rates through the whole 
periodic table, we have also calculated muon capture widths for a few selected nuclei, which will be also studied below 



It is approximately described by the Michel distribution 
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TABLE III: Experimental and theoretical total muon capture widths for different nuclei. Data are taken from Ref. :5f)|, and when 
more than one measurement is quoted in l5^ . we use a weighted average : r/o-^ = ""'^^^ ^/""^ = J2i ^/'^i- Theoretical results 

have been obtained by using non-relativistic kinematics for the nucleons (Sect.|n]of the Appendix). To illustrate the role played by the 
RPA correlations, we quote two different theoretical results: i) Pauli+Q obtained from Eq. 1621 without including FSI effects and RPA 
correlations (i.e., replacing j4^p^ by Af" in Eq. 1621 '). but taking into account the value of Q; ii) the full calculation, including all nuclear 
effects with the exception of FSI, presented in SectE] denoted as RPA. Finally, in the last column we show the relative discrepancies 
existing between the theoretical predictions given in the third column and data. 



|fc| [MeV] 

90 85 80 75 70 




FIG. 8: Inclusive muon capture differential width from ^^C, as a function of the outgoing neutrino energy (top axis) and of the energy 
transfer (bottom axis). Non-relativistic kinematics has been used for the nucleons. The two calculations are labeled as in Table ITTTl 

in Sect. IVIBl Our results are compiled in Table ITTll Data are quite accurate, with precisions smaller than 1%, quite 
far from the theoretical uncertainties of any existing model. Mediurn polarization effects (RPA correlations), once 
more, are essential to describe the data, as was already shown in Ref. Despite the huge range of variation of the 
capture widths^^, the agreement to data is quite good for all studied nuclei, with discrepancies of about 15% at most. 
It is precisely for ^^C, where we find the greatest discrepancy with experiment. Nevertheless, our model provides one 
of the best existing combined description of the inclusive muon capture in ^^C and the LSND measurement of the 
reaction ^^C {vfi,^~)X near threshold fl8|. 

Finally, in Fig.|Slwe show the outgoing energy distribution from muon capture in ^^C, which ranges from 70 to 



21 Note, r™P varies from about 4x10* s"! in to 1300 xlO* s"! in ^OSpb 
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FIG. 9: RPA and Coulomb (CB) corrections to electron and muon neutrino and antineutrino QE cross sections for different nuclei, as a 
function of the neutrino energy. A relativistic treatment of the nucleons is undertaken and FSI effects are not considered. J?rpa+CB is 
defined as (ctrpa+CB — fo)/o"0) where ctq does not include RPA and Coulomb corrections, while (trpa+cb includes these nuclear effects. 



90 MeV. The energy transferred to the daughter nucleus (^^B) ranges from to 20 MeV. We also show in the figure 
the medium polarization effect on the differential decay rate. As already mentioned, the shape of the curves in Fig.|S| 
will significantly change if a proper finite nuclei treatment is carried out, with the appearance of narrow peaks, but 
providing similar values for the integrated widths |4S|. 



B. Inclusive QE Neutrino and Antineutrino Reactions at Intermediate Energies 

In this subsection we will present results on muon and electron neutrino and antineutrino induced reactions in 
several nuclei for intermediate energies, where the predictions of the model developed in this work are reliable, not 
only for integrated cross sections, as in the previous subsection, but also for differential cross sections. We will present 
results for incoming neutrino energies within the interval 150-400 (250-500) MeV for electron (muon) species. The 
use of relativistic kinematics for the nucleons leads to moderate reductions of both neutrino and antineutrino cross 
sections, ranging these reductions in the interval 4-9%, at the intermediate energies considered in this work. Such 
corrections do not depend significantly on the considered nucleus. 

In FigO the effects of RPA and Coulomb corrections are studied as a function of the incoming neutrino/ antineutrino 
energy. These corrections are important (20-60%), both for neutrino and antineutrino reactions, in the whole range of 
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FIG. 10: Electron and muon neutrino and antineutrino inclusive QE cross sections and ratios for different nuclei, at intermediate energies. 
Results have been obtained with the full model without FSI, and using relativistic kinematics for the nucleons. For comparison we also 
show results obtained in the free space (low density limit, LDT, Eq. I^ ). 
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FIG. 11: Muon neutrino and antineutrino relativistic QE differential cross sections from different nuclei and several f^,!^^ energies. 
Results, denoted as 'Pauli+Q' or 'Pauli+Q' have been obtained in ^-^C and do not include RPA, FSI and Coulomb effects, while the rest 
of results have been obtained with the full model without FSI. 



considered energies. RPA correlations reduce the cross sections, and we see large effects, specially at lower energies. 
The RPA reductions become smaller as the energy increases. Nevertheless for the higher energies considered (500 
and 400 MeV for muon and electron neutrino reactions, respectively) we still find suppressions of about 20-30%. 
Coulomb distortion of the outgoing charged lepton enhances (reduces) the cross sections for neutrino (antineutrino) 
processes. Coulomb effects decrease with energy. For antineutrino reactions, the combined effect of RPA and Coulomb 
corrections have a moderated dependence on A and Z. Coulomb corrections reduce the outgoing positive charged 
lepton effective momentum inside of the nuclear medium. Thus, the phase space correction factor \K,' {r)\E[{r) /\k'\Ei 
is smaller than one and the cross section gets smaller. This effect, obviously grows with Z. On the other hand, the 
RPA suppression decreases when the lepton effective momentum increases and it grows with A. The combined effect 
explains the nuclear dependence found in the antineutrino plots. At the higher energy end the A— dependence becomes 
milder, since Coulomb distortion becomes less important. In the case of neutrinos, the increase of the cross section 
due to Coulomb cancels out partially with the RPA reduction. Finally, the existing differences between electron and 
muon neutrino/ antineutrino plots are due to the different momenta of an electron and a muon with the same energy. 

In Fig.^1 we show electron and muon neutrino and antineutrino inclusive QE cross sections and ratios for different 
nuclei, as a function of the incoming lepton energy. Results have been obtained with the full model presented in 
Sect lllll including all nuclear effects with the exception of FSI, and using relativistic kinematics for the nuclcons. 
Neutrino cross sections scale with N (number of neutrons) reasonably well, while there exist important departures 
from a Z (number of protons) scaling rule for antineutrino cross sections. These departures can be easily understood 
from the discussion of Fig|3 To better disentangle medium effects, the free space neutrino/ antineutrino nucleon cross 
section multiplied by the number of neutrons or protons is also depicted in the plots. 

In Fig. 111! we show muon neutrino and antineutrino inclusive QE differential cross sections as a function of the 
energy transfer, for different isoscalar nuclei and different incoming lepton energies. We see an approximate A— scaling 
and once more the important role played by the medium polarization effects. Similar results (not shown in the figure) 
are obtained from electron species. 

The double differential cross section da/dE'id\q \ for the muon neutrino reaction in calcium is shown in Fig. ^] In 
the top panel, we compare the lepton scattering angle distribution for three different values of the energy transfer. 
As usual in QE processes, the peaks of the distributions are placed in the vicinity of \q \ = y'2AI<f. In the bottom 
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FIG. 12: Muon neutrino differential cross sections in calcium as a function of the lepton scattering angle (top) and of the momentum 
transfer (bottom). The neutrino energy is 371.8 MeV. Top: Cross sections, without FSI and using relativistic kinematics for the nucleons, 
at different muon energies. Crosses have been obtained without RPA and Coulomb effects, while the curves have been obtained with the 
full model (up to FSI effects). Bottom: Cross sections, for a muon energy of 228.6 MeV, obtained by using relativistic (long dashed 
line, 'REL') and non-relativistic nucleon kinematics. In this latter case, we present results with (solid line , 'FSI') and without (short 
dashed line, 'NOREL') FSI effects. For the three cases, we also show the effect of taking into account RPA and Coulomb corrections 
(lower Unes at the peak). The areas (in units of IQ-*" cm^/MeV) below the curves are 3.50 (REL), 3.87 (NOREL) and 3.77 (FSI) when 
RPA and Coulomb corrections are not considered, and 3.13 (REL), 3.49 (NOREL) and 3.53 (FSI) when these nuclear effects are taken 
into account. 



panel, we show FSI effects on the differential cross section for one of the energies {Ej^ = 228.6 MeV) studied in the 
upper panel. We also show the effects of using relativistic kinematics for the nucleons. As anticipated, FSI provides a 
broadening and a significant reduction of the strength of the QE peak. Nevertheless the \q \ integrated cross section 
is only slightly modified (a reduction of about 2.5% when RPA corrections are not considered and only about 1% 
enhancement when they are included). 

In Fig. 1131 we plot double differential cross sections for fixed momentum transfer, as a function of the excitation 
energy. We show neutrino and antineutrino cross sections from ^^O. FSI effects are not considered in the top panel, 
and one finds the usual QE shape, with peaks placed, up to relativistic corrections, in the neig hborhood of q'^/2M. 
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FIG. 13: Ue and Ue diflferential cross sections in ^'^O as a function of the excitation energy, for fixed values of the momentum transfer 
and -E^— = 400 MeV. Top: Results obtained from the full relativistic model without FSI, with ('RPA') and without RPA and Coulomb 
corrections ('Pauli+Q(Q)'). Bottom: Results obtained by using relativistic (long dashed line, 'REL') and non-relativistic nucleon 
kinematics. In this latter case, we present results with (solid line, 'FSI') and without (short dashed line, 'NOREL') FSI effects. For the 
three cases, we also show the effect of taking into account RPA and Coulomb corrections (lower lines at the peak). The areas (in units of 
10~*° cm^/MeV) below the curves are 1.02 (REL), 1.13 (NOREL) and 1.01 (FSI) when RPA and Coulomb corrections are not considered, 
and 0.79 (REL), 0.90 (NOREL) and 0.85 (FSI) when these nuclear effects are taken into account. 



Once more, medium polarization effects are clearly visible. FSI corrections are studied in the bottom panel, and we 
find the expected broadening of the QE peak, but the integrated cross sections remain almost unaltered. 

Finally, in Table ITvl we compile muon and electron neutrino and antineutrino inclusive QE integrated cross sections 
from oxygen. We present results for relativistic and non-relativistic nucleon kinematics and in this latter case, we 
present results with and without FSI effects. Though FSI change importantly the shape of the differential cross 
sections, it plays a minor role when one considers total cross sections. When medium polarization effects are not 
considered, FSI provides significant reductions (13-29%) of the cross sections j^^. However, when RPA corrections 
are included the reductions becomes more moderate, always smaller than 7%, and even there exist some cases where 
FSI enhances the cross sections. This can be easily understood by looking at Fig. ^1 where we show the differential 
cross section as a function of the energy transfer for — 375 MeV. There, we see that FSI increases the cross 
section for high energy transfer. But for nuclear excitation energies higher than those around the QE peak, the RPA 
corrections are certainly less important than in the peak region. Hence, the RPA suppression of the FSI distribution 
is significantly smaller than the RPA reduction of the distribution determined by the ordinary Lindhard function. 
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[MeV] 


<7(i«0(i/,„M-^)) [10-*" cm^] 


a(i«0(^p,M+X)) [10-*° cm^] 


500 Pauli + Q(Q) 
RPA 


REL NOREL FSI 
460.0 497.0 431.6 
375.5 413.0 389.8 


REL NOREL FSI 
155.8 168.4 149.9 
113.4 126.8 129.7 


375 Pauli + Q(Q) 
RPA 


334.6 354.8 292.2 
243.1 263.9 243.9 


115.1 122.6 105.0 
79.8 87.9 87.5 


250 Pauli + Q(Q) 
RPA 


155.7 162.2 122.5 
94.9 101.9 93.6 


63.4 66.4 52.8 
38.8 42.1 40.3 


[MeV] 


a(^''0{ue,e-X)) [10-*° cm^] 

V ^ / y L J 


a(i«0(Pe,e+X)) [10-*° cm^] 


310 Pauli + Q(Q) 
RPA 


REL NOREL FSI 
281.4 297.4 240.6 
192.2 209.0 195.2 


REL NOREL FSI 
98.1 104.0 87.2 
65.9 72.4 73.0 


220 Pauli + Q(Q) 
RPA 


149.5 156.2 121.2 
90.1 97.3 92.8 


60.7 63.6 51.0 

36.8 40.0 40.2 


130 Pauli + Q(Q) 
RPA 


37.0 38.3 28.8 
20.6 22.3 23.3 


21.1 21.9 16.9 
10.9 11.9 12.8 



TABLE IV: Muon (top) and electron (bottom) neutrino (left) and antineutrino (right) inclusive QE integrated cross sections from 
oxygen. We present results for relativistic ('REL') and non-relativistic nuclcon kinematics. In this latter case, we present results with 
('FSI') and without ('NOREL') FSI effects. Results, denoted as 'RPA' and 'Pauli+Q(Q)' have been obtained with and without including 
RPA correlations and Coulomb corrections, respectively. 



VII. CONCLUSIONS 

The model presented in this paper, which is a natural extension of previous works 0, 0, 0, on electron, photon 
and pion dynamics in nuclei, should be able to describe inclusive QE neutrino and antineutrino nuclear reactions at 
intermediate energies of interest for future neutrino oscillation experiments. Even though the scarce existing data 
involve very low nuclear excitation energies, for which specific details of the nuclear structure might play a role, 
our model provides one of the best existing combined description of the inclusive muon capture in ^^C and of the 
measurements of the ^^C {v^,^~)X and ^^C {ve.,e~)X reactions near threshold. Inclusive muon capture from other 
nuclei is also successfully described by the model. 

The inclusion of RPA effects, in particular the nuclear renormalization of the axial current, turned out to be 
extremely important to obtain an acceptable description of data. This had been already pointed out in Refs. |l2lll6l 
EM I2H , and it is a distinctive feature of nuclear reactions at intermediate energies 0, 0, IS • C)n the other hand 
FSI effects, though produce significantly changes in the shape of differential cross sections, lead to minor changes 
for integrated cross sections, comparable to the theoretical uncertainties, once RPA corrections are also taken into 
account. 

The natural extension of this work is the study of higher transferred energies to the nucleus, also relevant for the 
analysis of future experiments aiming at determining the neutrino oscillation parameters with high precision. For 
those energies, the production of real pions and the excitation of the A(1232) or higher resonances will be contributions 
to the inclusive neutrino-nucleus cross section comparable to the QE one, or even larger. 
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FIG. 14: Muon neutrino QE differential cross sections in oxygen as a function of the energy transfer. The neutrino energy is 375 MeV. 
We show results for relativistic (long dashed line, 'REL') and non-relativistic nucleon kinematics. In this latter case, we present results 
with (solid line, 'FSI') and without (short dashed line, 'NOREL') FSI effects. We also show the effect of RPA and Coulomb corrections 
(lower lines at the peak). The integrated cross sections can be found in Table ITvl 



APPENDIX A: CC NUCLEON TENSOR 



1. Impulse Approximation 

Performing the traces in Eq. I|27() and taking into account that in Eq. (|28|1 both the particle and the hole nucleons 
are on the mass shell [p^ = (p + ^ ^'^"^ ^ 2p ■ q + = 0), one finds 

A^'^ip, q) = f {(p + qTf + ip + qYp^ + ^.g^^} + 2qHpvFYr [^d'"' ~ ' ^ P'<l'' +fP'' 

- I'^l'^i^ + Jp)} - IGF^t^vF^iq^q'' - q^gn + ^Gl \^2p^ p- + q^^ p- + p^^ q- + g^^- (t. _ 2M^) 

The above tensor admits a decomposition of the type 

and from Eq. we have 



1 _ M2 
4 g2 



a,iq') = 32iFrr~8{nvFYr^+8G\ 
asiq^) = IGGAiF^" + txvFY) 
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2. RPA Corrections 

Taking q in the z direction and after performing the RPA sum of Fig. we find, neglecting^^ corrections of order 
O {kFP^/]VP,kFp''^/M^, kpq^/M) 

iM^ = 8(F, ) jCp, (^^j + — ^-2 — {^^vF,) | — + 



^RPA .rp^^2/-. W)2p^_tM^'?V\ 9% ^^.,2 f 2p, + \q\ , I , g°(2p, + |g |) 
iM^ = 4(^^i ) |Cn^^^^ M^/-mJ^^^^2 ) l^f^^ ^ + 2^ + ^2 

- ^F,^,vF, —^^2G^i^C^— ,^ + ^^-C^-^^\-^^+2j^ (A5) 



toJ — g^ \ — q 



+ 2Gi{c^+^i^} (AT) 
^RPA _ .^y , „ ^v^^9V ^ |g|g(p) 



4A.2 - ^^GAFr + .vF/ ) - J (AS) 

with the polarization coefficients defined as 

^^(^)=|f-co/'Jc/^(g,MP' ^^^^^=|f-c/(g,M^*(g)P' ^^^^^^|f-c/feM^(g)P ^^'^^ 

In order to preserve a Lorentz structure of the type g^q", for the pseudoscalar-pseudoscalar and pseudoscalar-axial 
vector terms of the CC nucleon tensor, we have kept the RPA correction to the term „[il^i ( m^'^-g^ ^ ^) ^rpaj 
despite of behaving like {q"/\q |)^ w 0{q^/]VP). 

APPENDIX B: BASIC INTEGRALS 

In a non-symmetric nuclear medium, the relativistic Lindhard function is defined as 

The two contributions above correspond to the direct and the crossed ph excitation terms, respectively. For positive 
transferred energy only the direct term has imaginary part, which is given by 

ImC7;^(g, fc?;, A;^) = / d'pTR{q,p, k"p, k^) = -M^^^^^^-^<d{E^i. - + q'^mE^p - £l){E^p - £l) (B2) 
J 27r|g| 



Note that q'^/AI is of the order \q\^/M^ and as mentioned in Sect. ITTl Al we have considered fiy °^ order 0(0). 
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with 



£'„ = Max |m. - + I"' I 'Z' ' I . = yS^TftiFT, (B4) 

being Max(...) the maximum of the quantities included in the bracket. To perform the dPp integrations in Eqs. H28|) 
and H62|l is important to bear in mind that, though the LFG breaks down full Lorentz invariance, one still has 
rotational invariance, thus we find 

r«(g, k^, fc^) = J d?pTH{qJ, kl, ep)E{p) ^\{E'^F+ £l) X IniC7i?(g, k^j^Xp) (B5) 

ffl(g, kl, fcP ) - j d?pTR{qJ, kl, fc^)p = (^^luvVuiq. kl^Xp) + fcF, k^)^ q (B6) 

fc^) = y d^pTR{q,p, k^, kPp)E'{p) = 1 + (f P + SP^E'^) X ImI7H(q, fc]^, fc^ ) (B7) 

fc^^, fc^) = J d^pTniqJ, kl, kl)E{p)p= (^^0(9, k'^p^kl) + fc?;, fc^ 9 (B8) 

= 1 d'p:FR{q,pXF,k^F)py - ^^^Y^ + g'g^ = 1,2,3 (B9) 



with 

a,?(g,fc^,A:^) = fc?,. A:?,) - M^ImUflCq, fc]^. A:?,) (BIO) 

= ^{g4lna[7^(5,fc?;,A:^)+4(g0)2i?00(g,fc^,fcP) + 4gVr^(9,fcF,fc^^ (Bll) 

APPENDIX C: NON-RELATIVISTIC REDUCTION OF THE RESULTS OF APPENDIX B 

We take a non-relativistic reduction of the nuclcon dispersion relation 

m^M+^^m (ci) 

which implies, for consistency, that in the definition of the imaginary part of the Lindhard function and in all integrals 
given in the Appendix IbI the factors M/E{p) and M/E{p + q) should be approximated by one. Thus, we have^^ 

Uiak- kn 2 f Qik-F~\p\mP + q\-k'^F) I . ^ . (C2) 
which correspond to the direct and the crossed ph excitation terms, respectively^^. For positive values of we have 

imu{q, k-p, fcp ) - / d^pHqJ, kiXp) = -Af2 Q(g°)Q(,~g') e(^j^ -E^p + q°)eiE^ - sniE^ - £n (C3) 

J 2Tr\q\ 



with 

1 

4^2 



:F{q,p,k-p,kPp) = -—eiq")Siq'' + m-E{p + q)mkF-\pm\p + q\-k^F) (C4) 

£P = Max 1^?, -qO^M+^(^-^4r)} , E'f" = M + ^-^^ (C5) 



^ ' 2Af V I? I 



2 y ^ 2M 



We suppress the subindex R to distinguish the new expressions from the former ones. 

For symmetric nuclear matter pp = pn = p, the above expression coincides, up to a factor two due to isospin, with the definition of C/jv 



given in Eq. (2.9) of Ref. |35 
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To perform the integrations implicit in Eqs. 128|l and (|62|l we need 

r°(g, fc^, kPp) = J d?pT{q,p, k],,kPp)E{p) - ^ (^F + £") X ImUiq, fc^, fc^) (C6) 

f{q,klXF) ^ J d'pHQ,P,k^,kPp)p= (^-\ + ^y^Uiq,k"p,k''j,)q (C7) 

i?o°((?,fc^,A;^) = j d''pT{q,p,kl,kl)E\p) = i + X ImC7(g,fc^,A:^) (C8) 

Riq.k^F.kl) - 1 d=^p^(g,p,fc5^,A:?,)^(p>-'- - T^^iq^k^.k^) q (C9) 



fc]^, fc^) = / d^pn<l,P, k^,kPp)pY = + ^-r^ <?*g^ J = 1, 2, 3 (CIO) 



with 



a{q,k'l„kPp) = 2A^f{T°(g,fcJ,fc^)- AfImC/(q,fc^,A:^)} (Cll) 
4|g 



h{q,klXF) = -^{2Mq°-\q\^)\uvU{q,kl,ep) (C12) 



APPENDIX D: FREE NUCLEON CROSS SECTION 



The cross section for the process vi + n ^ I + p is given by 

'^-J = ^^TT: 77^ / dq^L^.A''^' 



87r(s - Af 2)2 



9mi„ 



p=(M,0) 



(Dl) 



where the leptonic {L) and nuclcon (A) tensors are defined in Eqs. and (|27IA1() . respectively, g^in(max) ~ ~ 

2Etj(Ei±\k' I) with E^, and , k' the incoming neutrino LAB energy and outgoing Icpton LAB energy and momentum, 
and finally s — {2E^ + M)M. The variable q^ is related to the outgoing lepton LAB polar angle {6') by q^ = {k — k')'^ = 
~ 2Ey{E[ — \k' I cos^^'). The tensor contraction in Eq. (|D1|I gives in the LAB frame: 



(g2 - r77|) \ ai + - - a4^| + (s - M2) j^-^aa - g'as }> (D2) 



p=(M,o) "12 2 2J ' ' [ 2 



p— (^iw.u; 

with the nucleon structure functions, ai{q'^), given in Eq. (|A3I) . 

The cross section for the process Di + p ^ + n is obtained from Eqs. ljDl|) and IID2|) by replacing 03 by —03. 
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